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Abstract 

2 ' In this paper we derive a bound on the rate of change of the gravita- 

tional constant G coming from the pulsating white dwarf G117-B15A. 
This star is a ZZ Ceti pulsator extensively studied with astroseismolog- 
ical techniques for last three decades. The most recent determination of 
P = (2.3 ± 1.4) x 10 -15 ss^ 1 for the 215.2 s fundamental mode agrees 
very well with predictions of the best fit theoretical model. The rate 
of change of the oscillation period can be explained by two effects: the 
cooling (dominant factor) and change of gravitational binding energy 
(residual gravitational contraction). Since the white dwarfs are pulsat- 
ing in g-modes whose frequencies are related to the Brunt- Vaisala fre- 
quency (explicitly dependent on G) observational determination of the 
change of the period (more precisely the difference between observed 
and calculated P) can be used to set the upper bound on the rate of 
change of G. In the light of the current data concerning G117-B15A we 
derive the following bound: |^| < 4.10 x 1CP 10 yr _1 . We also demon- 
strate that varying gravitational constant G does not modify cooling 
of white dwarfs in a significant way at least at the luminosities where 
white dwarfs are pulsationally unstable. 

Key words: gravitation - stars:oscillation - white dwarfs - stars:individual:G117- 
B15A 
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1 Introduction 



There is a renewed debate in the literature over the issue whether the quantities 
known as the constants of nature (such like G, c, h or e) can vary with time 
[T] . One of the reasons for this debate is associated with the advances of string 
theory and associated ideas that the world we live in may have more than four 
dimensions. 

String theory 2J is the only known framework which gives hopes to reconcile 
gravity theory with quantum mechanics and one of its hints is that the coupling 
constants appearing in low-energy Lagrangian are determined by the vacuum 
expectations of some a priori introduced massless scalar fields (dilaton or mod- 
uli fields). The interest in physical theories with extra spatial dimensions has 
also experienced considerable revival. In the framework of multidimensional 
theories [3] the four dimensional gravity constant, for example, is in fact an 
effective one appearing when integrated over additional dimensions. 

Another trigger for the interest in varying the fundamental constants is 
associated with recent observational evidence that the fine structure constant 
may indeed vary with time jl]. 

All this motivates us to take seriously the possibility that (at least some 
of) the fundamental constants could in fact vary in time. Although from his- 
torical perspective the gravity constant G was the first fundamental constant 
considered as a dynamical quantity 5J, in fact it is not strongly constrained 
concerning its variability nor even measured with accuracy compatible to other 
fundamental constants — latest CODATA report raised the relative uncer- 
tainty of G from 0.013% to 0.15% 0. In the past, there were numerous 
attempts to constrain the possible temporal variability of G as summarised 
in a recent review paper by Uzan 7j. White dwarf cooling has already been 
used to constrain the temporal variation of G in [§] at the level of white dwarf 
luminosity function. 

This paper is devoted to another constraint on ^ obtained from astero- 
seismological considerations of the ZZ Ceti pulsating star G117-B15A. More 
precisely, recent measurements of the rate of change of period of 215.2 s pulsa- 
tional mode [H] turned out to be in a very good agreement with evolutionary 
models of DAV white dwarfs comprising cooling and gravitational contraction 
effects. Pulsations of DAV white dwarfs are excited as the so called g-modes 
which arise due to the competition between buoyancy and gravity forces acting 
on matter element displaced in a stratified medium. We take an advantage 
of this circumstance to derive a bound on the rate of change of the gravity 
constant. 
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2 The role of changing G in white dwarf evo- 
lution 



White dwarfs are the final stages of the evolution of stars with initial masses 
of up to several M after they ascended the asymptotic giant branch, had 
lost their envelopes in the planetary nebula stage and developed a degenerate 
carbon-oxygen cores which is an internal energy reservoir and a thin non- 
degenerate envelope made up of He and H shells. Their subsequent evolution 
is governed by cooling, first dominated by neutrino losses throughout their vol- 
umes, later by surface photon emission and residual gravitational contraction. 
At certain stages of their evolution white dwarfs become pulsating stars (e.g. 
ZZ Ceti variables) and their pulsation period increases as the star cools down 

[101 mill- 
Under assumptions valid for DAV pulsating white dwarfs to which G117- 
B15A belongs, i.e. that neutrino emission is negligible and the star is in the 
evolutionary stage prior to cristallization of its core, the luminosity reads: 

M Ac r M A v A n ~l 

T ^dm = - (^ + p^_) rfm (l) 
o dt Jo dt dt ' y ' 

where s, u, P and p have usual thermodynamical meaning with lower case 
letters referring to specific (per unit mass) quantities. 

The second term in is in fact the work done on the star by gravitational 
compression. Our goal now is to see how the assumption that G is varying in 
time modifies the compression term. 

For this purpose let us start with hydrostatic equilibrium equation 

dP Gm 
dm ~ 47rr 4 

Multiplying (j2j) by 4irr 3 and integrating by parts one gets the virial theorem 
in a standard form: 

r M p 



Vl = -3 —dm (3) 
Jo p 

Differentiating gravitational energy Q with respect to time one obtains: 

rM p rM p 

Q = -3 / —dm + 3 —pdm (4) 
Jo p Jo p 2 

Then the term —3 J M — dm can be calculated by differentiating hydrostatic 



equilibrium equation (J2J) with respect to time, and multiplying it by Anr 

3 dP ^Gm . G Gm 
dm r 2 G r 
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On the other hand: 



fM Gm , G„ 
1 r z G 



Integrating (JSj) by parts over m and using (JEJ), one gets: 

-3 / - dm = AVt - 3— ft 

Jo p G 

Substituting this to (@J) one arrives at the formula: 

G f M P 
Q = —Q — / — p dm 
G Jo p 2 



This means that the compression term in $1} is equal to 



M do' 1 ■ G 

P ^L_ dm = n-^n (7) 

o at G 
Hence the luminosity of a white dwarf reads 

dE G ^ ... 

where: E = U + Q denotes total energy of the star and respective formula in 
Garcia-Berro et al. is recovered. 

However this does not mean that varying G can influence white dwarf 
cooling in a significant way (at least during most of their evolution). In order 
to see this, let us consider total internal energy of the star U as a function of 
state U (p, T) 

dU r M du , [ M f dudT dudp^ , 
dm = / (— — h — — — J dm 



dt Jo dt Jo dT dt dp dt 

using well known identity: — p 2 §^ = — P one gets: 

dU f M dudT j_ f M ^dPdp- 1 f M ^dp" 1 dE thermal , dE grav <Kl , G 



, dm+ T—^—dm-\ P^— dm = 

dt Jo dT dt Jo dT dt Jo dt dt dt dt G 

rearranging this we obtain the following: 

dE _ dC/ dn _ dE thermal dE grav G 

dt ~ dt dt ~ dt dt G 

After substituting this result to (JSJ) one can see that G terms will cancel leading 
to the standard formula: 
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dE- ^'thermal dEg rav 



dt dt 

Here, as commonly accepted, E therma i denotes thermal energy of the star, 
and E grav - the fraction of gravitational energy Q contributing to the luminos- 
ity. 

Although formula (jH)) is formally correct, yet the cancellation of G terms is 
still present in it in a hidden form. For this reason, basing on a standard lore, 
one is tempted (as we indeed were) to identify the change in the total energy 
of the star ^| with the change of its thermal energy dEth ^mai without noticing 

that the part of ^| coming from G gets cancelled with ^Vt term in 

Because dE ^ ma i = f Q M Cv <§ dm and = J M T§^^± dm one can 

see immediately that (even with changing G) zero temperature "black dwarf 
configuration would have zero luminosity as expected. 

In order to assess the real contribution of varying G to the white dwarf 
cooling let us assume that we have highly degenerate electron gas, ideal gas of 
ions and Coulomb correction: P = P e + Pj + P c , P e = ^ZifiPi, P c = — 0.3rPj 
where ip is degeneracy parameter and as usually T measures Coulomb potential 
energy compared with thermal energy. Then to the first order in 

dt Jo dt Jo Pi + Pc + Pe dt K J 

M P dp~ l , , , 

dm = (10) 



o 1 - 0.3r + OAZip dt 

= ((1 -o.3r + o.4ZV) _1 )(fi- S^) (ii) 

where the average above is defined by equation. Now, in carbon-oxygen white 
dwarfs like G117-B15A typically ip ~ 550 and V pa 55. Hence the modification 
to star's luminosity coming from varying G is about three orders of magnitude 
smaller than in the case of naive interpretation of formula (jSJ). At low lumi- 
nosities, however, close to the luminosity function cut-off the effect of varying 
G could be non-negligible as pointed out in 9J. However, in light of both the- 
oretical and observational uncertainties in precise determination of the cut-off 
the accuracy of luminosity function method is not a great one. 

Although changing G has no direct effect on cooling properties of white 
dwarfs during most of their evolution there hopefully exists certain possibility 
to obtain such a bound from DAV pulsating white dwarfs. The basic idea 
behind is simple: because DAV oscillations are driven by gravity force (acting 
against bouyancy force) the resulting oscillation period can be used to "mea- 
sure" G in a similar manner like in an elementary experiment one measures 
the acceleration of gravity g from measuring the period of a pendulum. 
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Extensive theoretical studies of non-radial oscillations of white dwarfs had 
been carried out long time ago [TT] . Appropriate theory of stellar oscillations 
consists in linearizing the Poisson equation as well as equations of momentum, 
energy and mass conservation with respect to small non-radial perturbations 

Hani. 

DAV pulsating white dwarfs oscillate in non-radial g-modes whose main 
restoring force is gravity (acting against buoyancy force). In spherically sym- 
metric stars, g-modes can be represented as standing waves /^ m (r, 9, 0, t) = 
3it!m( r )^i m (^ (P) ex P( — i&kimt) where by /' we denoted Eulerian perturbation 
of given quantity e.g. the pressure or density. G-modes observed in DAV pul- 
sators have modest radial orders 1 < k < 25 and low angular degrees 1 < I < 2 
(the latter is probably an observational selection effect). Because buoyancy is 
the restoring force for g-modes, the Brunt- Vaisala frequency N is the most 
important quantity setting the scale in the pulsation spectrum. Indeed in the 
asymptotic theory, eigen frequencies of g-modes are proportional to N. The 
Brunt-Vaisala frequency is usually expressed as: 

^_ 9 ^_ 9( ^_i-E) (12) 

where: g denotes local gravity, p - density, r - radial coordinate, p is the 
pressure and Ti is the adiabatic index |14j . In the context of white dwarfs 
alternative expression given by Brassard is even more useful 

iV 2 = — — (V ad -V + i?) (13) 
V Xp 

where X a := a = T, p and B = -± £^1 with Xl denoting the 

mass fraction of atomic species i. From the above formula one can see that 
the C/O, He/C and H/He transition zones manifest themselves as three peaks 
in N 2 which play an important role in modifying the period structure leading 
to the phenomena of mode trapping (near the surface) or mode confining (at 
central parts of the star) j27] . In fact similar asteroseismological considerations 
have already been used to calibrate compositional stratification of G117-B15A 
star |2Hl CHI 123- 

For a zero-temperature degenerate electron gas A = and hence (Jum = 
meaning that no g-modes are supported. However, if non-zero thermal effects 
are taken into account one can show [TTj that A oc T 2 and consequently o\ lm oc 
T 2 i.e. the periods scale like 1/T where T is the temperature of isothermal 
core. For high overtone modes one has: 

2 K_r_ 

Agl(l + 1) { ] 
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Following Winget, Hansen and Van Horn ^E] (see also [HI El) one can 
write: 

dlnP dlnT , dlnR , , 

-^ = -a^— + 6^— (15) 

dt dt dt y ' 

where: coefficients a and 6 are positive numbers of the order of unity. Since our 

further calculations will be an order-of-magnitude estimates we shall assume 

the a and h coefficients as equal to 1. First term comes from white dwarf 

cooling and dominates over the term from gravitational contraction. In fact 

DA white dwarf models (representative for ZZ Ceti pulsators such like G117- 

B15A) give a core temperature of 1.2 x 1§ 7 K and a cooling rate of 0. 005.fr 

i 

In the case of varying G formula ()16j) gets modified 

dlnP dlnT ^^dlnR dlnG , 

dt dt dt dt 

This means that time varying G affects the contraction term which is in agree- 
ment with underlying physics as described above. 



3 Asteroseismology of G117-B15A and bound 
on variability of G 

Since its discovery in 1976 (201 G117-B15A has been extensively studied and its 
basic physical properties are quite well established. Regarding its variability 
the observed periods are 215.2, 271 and 304.4 s together with higher harmonics 
and linear combinations thereof [21]. Quite recently, with a considerable time 
interval of acquired data, Kepler et al. [S] recalculated the rate of period 
increase and found a value of P = (2.3 ±1.4) x 10~ 15 s s^ 1 . It has been 
claimed JH] that the 215.2 s mode of G117-B15A is the most stable oscillation 
ever recorded in the optical band (with a stability compared to milisecond 
pulsars) and that observed P value is really due to evolutionary effects. This 
circumstance has already made it possible to derive bounds on axion mass ^H] 
and the compactification scale M s in the Arkani-Hamed-Dimopuolos-Dvali 
theory with large extra dimensions [2*2] . 

The white dwarf pulsator G117-B15A has a mass of 0.59 M (established 
spectroscopically) , effective temperature T e // = 11620 K [2Sj and luminosity 
log{C/C Q ) = -2.8 [21 (i.e. £ = 6.18 x 10 30 erg s" 1 ). More recent study 
by Koester and Allard [22] suggested the lower value for the mass 0.53 M Q . 
Typical model for such CO star (with C:0 = 20:80 [27J or C:0=17:83 [2H]) 
predicts the central temperature T = 1.2 x 10 7 K JH] and the radius R = 
9.6 x 10 8 cm 8j. Knowledge of the internal chemical profile of the core 
is essential for the following reasons. First, the heat capacity of the core is 
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mainly due non-degenerate ions and is inversely proportional to the average 
atomic number, thus enhanced O abundance implies less heat capacity, faster 
cooling and hence larger P. Second, the smaller C abundance implies steeper 
abundance slopes at the outer boundary of the core. This can lead to different 
radial configurations of the modes. It turns out that the 215.2 s mode has 
large amplitude in this interphase and can be used to constrain the admissible 
range of C abundance |18j . The chemical composition of the white dwarf 
interior as a function of mass calculated by Salaris et al. [2E] are in agreement 
with Bradley's results concerning DA pulsators j2Zj as well as with the work 
of Dominguez et al. who computed updated chemical profiles for the AGB 
progenitors of white dwarfs finding similar stratifications [26J. Therefore they 
can serve as a reliable reference. 

There are two main processes which govern the rate of period change in 
theoretical models of ZZ Ceti stars: the cooling of the star (which increases the 
period as a result of increasing degeneracy) and residual gravitational contrac- 
tion (which shortens the period) [TT]. For the G117-B15A star the contraction 
rate is negligibly small [B] hence confirming the assumption that cooling is 
a dominant process for this star. Theoretical value of P used below in our 
estimates was derived from numerical calculations JH] m which such issues 
as gravitational contraction, chemical composition of the star or properties 
of outer He-H layers were carefully taken into account. This value is also in 
agreement with previous independent evolutionary calculations 29 . There- 
fore, the correction from residual gravitational contraction is already present 
in theoretical rate of period change adopted here. 

The difference between observed (labeled " O" ) and calculated (labeled " C" ) 
logarithmic time derivatives of the period (precisely an absolute value thereof) 
|A(£)| := |(|)o - (|)o| reads: |A(|)| = |g|. So the upper bound for |A(|)| 
translates into a bound on the rate of change of G. 

The most recent determination |8] of P for the 215.2 s mode gives the 
value (2.3 ±1.4) x 10~ 15 s s~ l whereas theoretical prediction for the rate 
of change of the period is 3.9 x 10~ 15 s s~ l 18J. A conservative method of 
constraining the admissible role of varying G in the discrepancy between theory 
and observations comes from the fact that the theoretical value falls within a 
two-sigma interval for the observed value. Therefore current observational 
knowledge concerning G117-B15A pulsating star and associated consequences 
for the P of the 215.2 s mode can be translated into the following bound for 
the rate of change of the gravitational constant: 

G 

| — | < 4.10 x 10" 10 yr" 1 (17) 

Gr 

Of course, while constraining the variability of G in a manner described 
above we have neglected other non-standard effects. For instance a larger 
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(than currently thought admissible) ammount of axion JH] o r Kaluza-Klein 
graviton [22] emission could in principle be compensated by respective varia- 
tion of G. Not having in mind any specific theory of varying G (nor any other 
exotic physical theory) we disregarded the issue of joint effect of varying G 
and extra cooling by emission of exotic particles. Of other effects influencing 
the mode properties in a way that may interfere with changing G, the begin- 
ning of cristallization would slow down the cooling rate allowing for larger G. 
However, it has been argued jHO] that G117-B15A is not cool enough to have 
a crystallized core. 

4 Conclusions 

It would be instructive to compare the derived bound on the rate of change of 
the gravitational constant with already existing bounds. They can be roughly 
divided into four classes: geophysical bounds, bounds from celestial mechanics, 
astrophysical and cosmological bounds. The extensive discussion of these issues 
can be found in a review paper by Uzan 7j. Let us here invoke only the 
strongest (or the most representative) bounds from these four categories and 
discuss their relation to our bound. 

In the class of geophysical/paleonotological bounds the strongest one is 
the bound -g < 8 x KT 12 yr 1 coming from the stability of radii of the 
Earth, the Moon and Mars as estimated by McElhinny |31| . Celestial me- 
chanics arguments in the Solar System refer either to systematic deviations 
from Keplerian orbital periods in the case of time varying G (which can reveal 
themselves e.g. in lunar occultations) or to monitoring of the separation of 
orbiting bodies (lunar laser ranging or planetary ranging). The most stringent 
bounds obtained by these techniques are: g < (3.2 ± 1.1) x 1CT 11 yr' 1 [32] 

and % < (-2. ± 10.) x HT 12 yr' 1 [33], respectively. 

In the case of pulsar timing data Kaspi et al. [31] derived the strongest 
bounds in this class for PSR B1913+16 as: § < (4. ± 5.) x 1(T 12 yr' 1 . 

Cosmological bounds come out of the big bang nucleosynthesis (BBN) con- 
siderations and the most stringent ones are those of Accetta et al. [HHI (|§| < 
9. x 10~ 13 yr- 1 ) and Rothman and Matzner EH (||| < 1.7 x 1(T 13 yr- 1 ). 
They are however valid only for the Brans-Dicke theory. Moreover, the BBN 
arguments based on the abundance of light elements are sensitive not only to 
the gravity constant but also (and to greater extent) to fine structure constant 
and other factors like the neutron to proton ratio, the neutron-proton mass 
difference or the number of neutrino families. 

The most stringent astrophysical constraint |g| < 1.6 x 10~ 12 yr- 1 comes 
from helioseismology [3*3] . However it has been derived for a specific Brans- 
Dicke type theory with varying G. Previous white dwarf constraint based on 
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a position of the cut-off on the luminosity function of white dwarfs |9J was 
-g < 3.±\ x 10- 11 yr-\ 

In conclusion one can say that the bound derived here for G from the 
observed rate of change of the pulsational period of G117-B15A supports the 
claims first expressed in JH] that this star is becoming a tool for fundamental 
physics. Moreover, there are other DAV pulsating stars such like L19-2 or 
R 548 jHH| for which the rate of change of the period can be measured in 
forthcoming WET campaigns. Having high precision asteroseismological data 
for a larger set of objects one would be able to test exotic physical theories 
(such like those in which G varies) with much better accuracy. 
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